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Abstract. Measurements of force generated by a positive unipolar wire-to-plane corona discharge in air 
are compared with numerical simulations. The generated force does not depend on the ion or electron mo- 
bilities, preventing the influence of uncertainty and variation of these parameters. A method is described 
to simulate the voltage and charge distribution for a wire-to-plane set-up. This method enables the deter- 
mination of the transition between unipolar and bipolar discharges. In the experimental set-up breakdown 
electric field of air reduces with increasing discharge current. 



PACS. 02.60.Cb Numerical simulation; solution of equations 
52. 80. He Glow; corona 



02.70.Bf Finite-difference method 



1 Introduction 

Atmospheric pressure corona discharges in air are impor- 
tant for use as electrostatic precipitators, ozone genera- 
tors [1], thunderclouds [2], mass transport [3] and so on. 
An upcoming field is the use of corona discharge to gener- 
ate electric winds for cooling purposes [4], The physics of 
corona discharge however is not fully understood. Peek [5] 
was one of the first to study the phenomenon in a system- 
atic way. He found the onset breakdown electric field at 
the wire surface of a wire-cylinder discharge corresponds 
to 

E a = fE s(l + ^\ , (1) 

where / is an irregularity factor (for polished copper wire 
/ w 1), E the breakdown electric field of air equal to 31 
kV/cm, S the relative density of the air compared to the 
density of air at 25 °C and 1013 kPa, a the wire radius 
and a a constant equal to 0.0949 cm. In this paper the 
outer edge of the corona is defined as the distance to the 
center of the wire, r c where the electric field equals the 
breakdown electric field. In the corona between wire sur- 
face and outer corona edge it is assumed that space charge 
has negligible influence on the electric field. This is a valid 
assumption as close to the wire the electric field gradient 
is already very large due to geometry. Then, according to 
equation (1) the corona size is given by 



(2) 



This can be calculated by using equation (11) with a(x) = 
1. The second term at the right hand side of this equa- 
tion represents the ionisation sheath. In this sheath the 



discharge is carried by electrons, positive and negative 
ions [6] . Outside this sheath the discharge is mainly unipo- 
lar. A recent overview of possible electrical breakdown cri- 
teria has been given by Lowke [7]. He couples a a to the 
discharge parameters of air 



HQ) 



(E /N ) 2 N B 



(3) 



where Q is the number of ionisations to occur to sustain 
breakdown, N a is the air number density at 25 °C and 
1013 kPa and B is 1.35 XlO 16 \- 2 mT 2 \ Q depends on 
the corona discharge mechanism. He suggests a possible 
mechanism for the onset corona discharge but does not 
take into account the effect of space charge in the air sur- 
rounding the corona. Normally this is not a problem as 
measurements of the breakdown electric field of air are 
mostly done for relatively small space charges at rela- 
tively small discharge currents. However, for the above- 
mentioned applications the breakdown electric field in the 
presence of a discharge current is important. Additional 
problems are the relative uncertainties of the mobility of 
the involved ions and electrons. In the following it is shown 
that the force generated by a wire-to-plane corona dis- 
charge is independent of charge carrier mobilities. This 
force only depends on the applied potential difference be- 
tween wire and plane and the geometry of the discharge. 
Force measurements were performed to determine the in- 
fluence of space charge on the breakdown electric field of 
air at atmospheric pressure and to eliminate the influ- 
ence of the charge carrier mobilities. In the following this 



1 Lowke uses B=2.08x 10 16 V^m -2 , but he takes 25 kV/cm 
for the breakdown electric field of air. 
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method is described and measurements for several geome- 
tries are presented and discussed. 



2 Outline of the method 

The method is based on solving the Poisson equation for 
certain geometry combined with the continuity equations 
for ion and electron flow. Solving Poisson equation gives 
the potential distribution, U and the continuity equations 
give the space charge distribution, p. In air (the relative 
permittivity, e r is approximated by 1) the Poisson equa- 
tion is 

V 2 [/ = — P - , (4) 

where e Q is the permittivity of free space. For each charge 
species i, being positive or negative ions or electrons, the 
stationary continuity equation is 



V.Ji = 



(5) 



where Ji is the current density of charge species i. It is 
assumed that no transformation from one species into an- 
other occurs. If the mobility, pi of the species does not de- 
pend on either space charge nor electric field, E = — V£7 
it is given by 

Ji = Pi\pi\E: . (6) 

where pi is the charge density of species i. Hence, equa- 
tion (5) reduces to 

PiV.E-1-E.Vpi = . (7) 

Adding all equations (7) results in the continuity equation 
for the total space charge 



pV.E + E.Vp = 



(8) 



The force exerted by the electric field on a point particle 
with charge q is by definition P = qE. This can be ex- 
tended to find the force exerted by the electric field on a 
charge density distribution in a volume V: 



pEdV 



(9) 



Neither equation (4) nor (8) contains any reference to 
the mobility of the involved ions or electrons, other then 
that they are constant. The simultaneous solution of these 
equations results in voltage and charge density distribu- 
tions independent of these parameters. Hence according 
to equation (9) the force exerted by the charge distribu- 
tion does not depend on the mobilities. Further, this force 
is uniquely determined by the applied voltage difference 
between wire and plane and the geometry. The electric 
current, I through some flat area A is given by 



Mi I A 



lE.n dA 



(10) 



where n is the normal to area A. The current consists of 
contributions from all charge species. Hence, it is much 
more difficult to calculate accurately than the force. Fur- 
ther, if the total charge distribution consists of positive 
and negative species, the current increases as indicated 
by the absolute sign. 
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Fig. 1. Sketch of a cross section of the wire-to-plane geometry 



3 Experimental set-up 



The detection of force exerted by a space charge distri- 
bution needs equation (9) to produce a result different 
from 0. This can be achieved by an asymmetrical corona 
discharge, like a wire-to-plane discharge. Normal the wire 
is put above and parallel to the plane. Then, the electric 
wind generated by the space charge will push against the 
plane and reduce the measurable force [8]. To overcome 
this problem, the plane is reduced to a sheet and the wire 
is put parallel to the sheet at a distance S. The wire is po- 
sitioned in the same plane as the sheet. H is the height of 
the sheet (4.0 cm in the experiments). This configuration 
is shown in figure 1. Sheet thickness was 2 mm, far less 
than either S or H. The upper and lower end of the sheet 
were smoothly rounded with a radius of 1 mm to as much 
as possible avoid production of opposite signed charges. 
The length of both wire and sheet was 53.5 cm. In the ex- 
periment S was varied between 2.0 and 6.0 cm with steps 
of 1.0 cm. Wires of different radius a were used: 0.040, 
0.0635, 0.075, 0.10, 0.1575 and 0.25 mm. Measurements 
were performed at atmospheric conditions: 20 ± 1 °C and 
1013 ± 10 kPa. Relative humidity was between 0.5 and 
0.8. The force was measured by the reduction of weight 
of the set-up with a lever. One arm of the lever was con- 
nected with light non-conducting wires attached to the 
set-up. The length of these wires was about 50 cm ensur- 
ing the influence of the surroundings was negligible. The 
other arm was rigidly connected to a balance with a range 
of 100 gram and an accuracy of 0.1 gram. The ratio of 
the length of the arms of the lever was 1:10. This lever 
enabled force measurements between and 100 mN with 
an accuracy of 0.1 mN. The measurements were done with 
increasing and decreasing voltages to check for systematic 
deviations. A DC high voltage between and 30 kV was 
applied to the wire while the sheet was grounded. An am- 
pere meter monitored the current supplied by the high 
voltage power supply. 
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4 Simulation 



To simulate the electric field a model was developed fol- 
lowing [9]. First, the problem is reduced to 2 dimensions 
by assuming an infinitely long wire and sheet. Second, 
the problem is made periodic in both x- and ^-directions 
with a period of 2W X and W y respectively. 2W X respec- 
tively W y being the distance in x-respectively y-direction 
between two adjacent set-ups. x — W x and the x = 
are symmetry axis. If W x and W y are large enough the 
results can be applied to a single set-up. Here, W x and 
W y are chosen optimizing the computation time and re- 
quired accuracy. Third, the problem is discretized by as- 
suming a square grid with spacing A. The number of grid 
lines in the x-direction, n equals W x /A + l, the number in 
the y-direction, m equals W y /A. Fourth, the appropriate 
boundary conditions arc defined and the discretized ver- 
sion of equations (4) and (8) are solved sequentially until a 
stable solution is found. For the solution of the discretized 
Poisson equation a Poisson solver was used based on the 
linear conjugate gradients method [10]. For the solution of 
the descritized continuity equation a non-linear conjugate 
gradient method was used [11]. The appropriate bound- 
ary conditions are not evident. At the sheet the potential 
was simply set to 0. At the wire position however, the 
wire radius had to be taken into account. As the wire ra- 
dius is very small compared with the grid size a method 
had to be devised to take it into account. The method 
as described in [9] was used and extended to incorporate 
unipolar space charge around the wire. Around the wire 
up to the first adjacent grid point cylindrical geometry is 
assumed. This is justified as A <C S. The space charge 
density at r = r c is taken to be p c . Then, the electric field 
and charge density near the wire are given by 
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E(r) = fE S 



a{r/r c ) 
r/r c 



p{r) = 



Pc 



a{r/r c ) 



(11) 



(12) 



where a(x) = \J\ + i](x 2 — 1) if x > 1 and a(x) = 1 if 
x < 1 and r\ — p c r c / fE Se is a measure for the total 
space charge at the outside of the corona (r = r c ). For 
a unipolar discharge (outside of the corona) r\ is propor- 
tional to the discharge current, / 



7] = 



I/L 



2^pe {fE 5f 



(13) 



where p is the mobility of the unipolar charge carriers. 
Note that r\ is determined by the experimental conditions 
only. A graph of the electric field as function of r for several 
77 is shown in figure 2. The theoretical maximum value for 
77 is 1. For 77 = 1 the electric field in the whole gap equals 
the breakdown electric field and there is no way to prevent 
spark discharge. For r/r c < 1 the lines coincide. Inside the 
corona the electric field gradient due to the space charge 
density is negligible compared to the geometric gradient. 
For all experiments 77 < 0.05, so the influence of the space 
charge on the electric field is negligible near the wire up 
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Fig. 2. Ratio of electric field, E(r) to electric field at corona 
edge, E(r c ) as function of distance to center of wire, r for 
77 = 0; 0.01; 0.1 and 0.9. The arrow indicates increasing values 
of Tj. 



to r = 5r c . Hence, the exact shape of the streamers in the 
corona area is not important for the results presented here. 
The voltage drop from wire to outer edge of the corona is 
given by 

AU c = -fE 5r c ln^ , (14) 

independent of r\. The voltage drop from outer corona edge 
to a radius of r > r c is given by 

AU(r) = -fE Sr c (a(-)-l + x(-)VT^) , (15) 



where 



X(x)=ln{ 

Vi - v + a \ x ) 



(16) 



This function is shown in figure 3 for different values of 77. 
Note the positive values for the voltage drop as r/r c < 1. 
Here, AU(r) = —fE 5r c ln(r/r c ), again independent of r\. 
For the boundary conditions at the wire the electric field 
at all grid positions adjacent to the wire is made equal to 
E(A) and the space charge density is set to p(A). 



5 Results and discussion 

The wire voltage was calculated by adding AU (A) + AU C 
to the average of the voltage at the wire adjacent grid 
points. From the simulations the spread in the voltages at 
these grid point was only a few percent. This indicates the 
accuracy of the assumption of cylindrical geometry near 
the wire up to the first grid point. / was determined by 
fitting the onset breakdown electric field to equation (1). 
For the onset r c was calculated with equation (2). As an 
example, the onset voltage distribution for S = 3.0 cm 
and a = 0.040 mm is shown in figure 4. The position of 
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0.88(2) 



Table 1. Fitted irregularity factor / and simulation parameters for all geometries. / was fitted to match the measured onset 
electric field at the wire position with equation 1. 





U I V 



Fig. 4. Voltage distribution, U at the onset voltage for S = 
3.0 cm, LS — 4.0 cm and a = 0.040 mm. Grid size was 5.0 mm; 
W x = 15.0 cm; W y = 30.0 cm. Wire onset voltage = 5.82 kV. 



the wire and sheet are evident. Notice the steep rice at 
the wire position. The wire onset voltage was 5.82 kV. 
In this case the corona radius according to equation (2) 
is 0.235 mm. The corona voltage drop according to equa- 
tion (14) was (/ = 1) 1290 V. In the simulation the voltage 
at the wire position equals the voltage at the outside of 
the corona. Hence, in the simulation the voltage at the 
wire position was 4.53 kV. The fitted irregularity factors 
are shown in table 1. The onset voltages varied between 
5 and 12 kV. The average irregularity factor is 0.86, this 
indicates that the used value of E Q is slightly too high. 
Using a value of 27 kV/cm instead of 31 kV/cm, increases 
the average irregularity factor to 0.99. The small spread 
in the fitted value of the irregularity factor is evidence for 
the correctness of the performed simulation. Using these 
irregularity factors, at every measured voltage and geom- 
etry the simulated force was fitted to the measured one 
by adjusting r c and p c . Both parameters are needed to 
determine the force at a certain voltage and geometry. All 
other parameters were taken from the experimental condi- 
tions. Figure 5 shows an example of a voltage distribution 
for p c = 128 Cm~ 3 (which corresponds to a unipolar 
discharge current of 0.1 mA/m, if fi = 2.0 cm 2 V _1 s _1 ) 
with the same geometry as in figure 4. Figure 6 shows an 
example of the corresponding space charge density distri- 
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Fig. 6. Charge density distribution, p for p c = 128 fiCm 3 
and same dimensions as in figure 4. Wire voltage is 13.6 kV. 
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Fig. 8. Discharge current, 7 as function of wire voltage, U for 
same dimensions as figure 4. The symbols with error bars repre- 
sent the measurements for increasing (crosses) and decreasing 
(circles) voltage. The lines represent the simulated results. 




5 10 15 20 25 30 

U I kV 

Fig. 7. Exerted force, F as function of wire voltage, U for same 
dimensions as figure 4. The symbols with error bars represent 
the measurements for increasing (crosses) and decreasing (cir- 
cles) voltage. The lines represent the simulated results. 



bution. Notice the small oscillations in the current den- 
sity. This is due to the limited accuracy of the calcula- 
tions and the accuracy of discretization. It is assumed 
they do not significantly change the results presented here. 
Simulations and measurements of the exerted force for the 
above-mentioned geometry as a function of applied wire 
voltage are shown in figure 7. The symbols with error bars 
represent the measurements for increasing and decreasing 
voltage. The lines represent the simulated results. The fit- 
ted force matches the measurement accurately. The dis- 
charge current is calculated using the voltage and charge 
distributions and equation (10) applied for a unipolar cur- 
rent with a mobility of 2.0 cm 2 V _1 s _1 . This is show in 



figure 8. For small voltages the deviation with the mea- 
surements are within the accuracy of the measurements. 
Note that both force and discharge current measurements 
are accurately calculated with the same space charge den- 
sity distribution. For the force this is a result of the fit. For 
the discharge current this is not evident, as the currents 
were not used in the fits. Force is a volume and current an 
area integral over a charge and electric field distribution 
corresponding to equation 9 and 10 respectively. Hence, 
good agreement between the measured and simulated dis- 
charge current is a further indication of the validity of this 
method. For larger voltage the measured discharge cur- 
rent grows faster than the simulated one. Here, the cur- 
rent measurements with increasing and decreasing volt- 
age starts to deviate slightly. This is caused by a change 
of the discharge from unipolar to bipolar. The charges of 
the other sign reduce the force but enhance the current. 
From the fitted r c the corresponding values of ^faZ are 
calculated using equation (2). They are shown in figures 9 
and 10. The vertical lines in the graphs show the max- 
imum force for which the simulated current matched the 
measured one. This is the maximum force for a unipolar 
discharge in the specified geometry. ^faZ decreases from 
the onset value of 0.0308 m 0,5 more or less linear with 
a slope dependent on the geometry. The slope is propor- 
tional to S. For different wire diameters the slope is ap- 
proximately the same. The reduction of ^faZ corresponds 
to a reduction of the electric field at the wire position. In 
literature this is attributed to a reduction of the break- 
down electric field of air dependent on the applied over 
voltage [12,13]. The physics behind this are not yet un- 
derstood. The reduction of yJaZ to corresponds, accord- 
ing to equation (3), to a reduction of Q to 1. The number 
of needed ionisations reduces to 1. The slope's indepen- 
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Fig. 9. Calculated ^faZ as function of measured force F for 
measured geometries with S = 2.0; 4.0 and 6.0 cm. Closed di- 
amonds: a=0.040 mm; closed squares: a=0.0635 mm; stars: 
a=0.075 mm; open diamonds: o=0.10 mm; open squares: 
a=0.1575 mm and dots: a=0.25 mm. Error bars are estimates 
of accuracy of simulation. The vertical lines indicate the max- 
imum force for which the discharge is unipolar. 
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Fig. 10. Calculated yfaZ as function of measured force F for 
measured geometries with S = 3.0 and 5.0 cm. See figure 9 for 
details. 



dence of a and its proportionality to S suggest Q depends 
on the total amount of charge present in the discharge 
and not on local conditions around the wire. Lowke pro- 
poses the breakdown of air is caused by ions produced by 
interaction of metastablc molecules. The effect presented 
here supports his conclusion, as the number of metastable 
molecules present in the air depends on the charge density. 
The measurements presented here arc not accurate enough 
to determine the mechanism of this dependence. For larger 
forces, where the discharge becomes bipolar, ^/aZ becomes 
negative, indicating the limit of equation (8). This equa- 



tion assumes no transformation of one charge species into 
another. This is further supported by the observed trend 
that for larger interaction times of the charge species (i.e. 
for large traveling distance of the charges) ^faZ becomes 
negative for smaller forces. 



6 Conclusions 

Force generated by an asymmetrical wire-to-plane dis- 
charge can be measured and accurately simulated. Com- 
parison of measurements of both force- voltage and current- 
voltage characteristics with finite differences simulations 
show good agreement for the unipolar region. These mea- 
surements can be used to determine the transition voltage 
of an unipolar into a bipolar discharge. Further, it is possi- 
ble to determine the breakdown electric field of air depen- 
dent on the discharge conditions. The presented measure- 
ments are in qualitative agreement with the breakdown 
mechanism as proposed by Lowke, but measurements are 
not accurate enough to determine the dependence of this 
mechanism on the space charge density distribution. 

The author wishes to thank Mrs. G. Gravendeel-van Ingen for 
a valuable contribution to the experimental set-up. 
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